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Abstract 

It is pointed out that the diagonal amplitude of the El — El chromo-electric inter- 
action with soft gluon fields (chromo-polarizability) can be measured directly for the 
J ftp and T resonances in the decays J /ip — > mrl + £~ and T — ► 7nrl + £~ with soft pions. 
For the J/V> this amplitude is often discussed in connection with the J/ip interaction 
with nuclear matter, while for the T the chromo-polarizability enters the estimates of 
the non-perturbative mass shift of the resonance relevant to precision determination of 
the b quark mass from the T mass. 



It is realistic to expect that a very high statistics data on heavy quarkonia will become 
available in a foreseeable future. In particular the plans of the CLEO-c experiment include 
acquiring about 10 9 events at the J/ip resonance and, given the capabilities of the B factories, 
a great increase in the statistics of the T resonances is also quite possible. With such amount 
of data quite rare processes can be studied, which can help in resolving some of the old- 
standing problems in dynamics of heavy quarkonia. The purpose of the present letter is to 
point out that a study of the decays J/ip — > ttit£ + £~ and T — > tttt£ + £~ with soft pions would 
allow to measure the gluon polarizability of the respective quarkonium resonance, which 
determines the strength of the quarkonium interaction with soft gluon field, and which thus 
far is only guessed theoretically. 

The chromo-polarizability a of a state of quarkonium can be defined (in complete analogy 
with the usual polarizability of atoms) through the effective Hamiltonian of the quarkonium 
state interaction with a soft chromo-electric field E a : 

H eff = - l -aE a -E a . (1) 

This effective interaction arises in the second order in the leading El chromo-electric dipole 
term in the multipole expansion [2J |3| of the interaction of a heavy quarkonium with the 
gluon field: 

H E1 = -^ a f-E a (0) , (2) 

where £ a = f{ — 1\ is the difference of the color generators acting on the quark and antiquark 
(e.g. t\ = X a /2 with A a being the Gell-Mann matrices), and r is the vector for relative 
position of the quark and the antiquark. (In the normalization used throughout this paper 
the QCD coupling g is included in the definition of the field, so that e.g. the gluon field 
Lagrangean reads as L = — (F^ v ) 2 /(kg 2 ).) Thus for a colorless S wave state the chromo- 
polarizability a is given by the diagonal matrix element: 

a = ^{S\CnG A nC\S) , (3) 

where Ga is the Green's function for a heavy quark pair in color-octet (adjoint) state. A 
theoretical understanding of such matrix element is at least highly model dependent at 
present due to its sensitivity to the wave function of the quarkonium and due to presently 
unknown propagator Ga of a colored quark pair. 

The chromo-polarizability of J/ip: aj/^, enters as an important parameter in analyses of 
the J/ip interaction with nuclear matter (see e.g. in Refs. [HIS]) and is estimated essentially 
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on dimensional grounds[6 . The same quantity for the T resonance: atr, determines [3] the 
non-perturbative shift of the T mass due to the gluon vacuum condensate[7] (0|(F^) 2 |0): 

5M x = ia T (0|(F;j 2 |0) . (4) 

The knowledge of the non-perturbative shift of the T mass is important for the method 
of determining the b quark mass by comparing My with the perturbative QCD expression 
in terms of m;, for the mass of the lowest 3 S\ bound state of the quark- ant iquark pair (for 
a detailed discussion see the review jH] and references therein). Practical estimates of this 
non-perturbative shift use an extrapolation down to the bottomonium of the formulas [§| ITU] 
for asymptotically heavy quarkonium. In the latter limit the heavy quarkonium is essentially 
a Coulomb system, and the matrix elements of the type in eq. (j3J) can be found explicitly. In 
particular for the IS* state the chromo-polarizability is found [HI US] as 

78 m Q 

ai5 = 425 ^' (5) 
with ks = 2rriQa s (kB)/3 being the Bohr momentum for the heavy quarkonium. If applied to 
bottomonium, this formula gives ay ~ lGeV~ 3 . Using then the value of the gluon vacuum 
condensates (0|(F^) 2 |0) = (0|4vra s (G^) 2 |0) « 0.5GeV A , one estimates from eq.® the 
non-perturbative shift of the mass of T: 5My « 50 — 60 Mel 7 . However, in addition to the 
uncertainty in the value of the gluon condensate, there are at least two other uncertainties 
involved in this estimate: one arising from the application of the asymptotic expression JSJ) 
to bottomonium, and the other associated with possible contribution of vacuum averages of 
higher dimension, generally nonlocal, gluonic operators. Although it is not clear at present 
to what extent the latter uncertainty can be estimated, a direct measurement of the chromo- 
polarizability ax would definitely fix at least one factor in this problem. 

It should be noted that the non-diagonal amplitudes of the type in eq. (J3J) for the transi- 
tions between the 25" and IS* states, 

ais-25 = ^(iS\CnG A nC\2S) , (6) 

can be found both in charmonium and bottomonium from the spectra and the rates of the 
pionic transitions ip' — ► TrnJ/tp and T' — ► ttti T. The relation arises through the fact that 
the dominant part of the amplitude of production of two pions by the gluonic operator (E a ) 2 
is determined [11 J by the trace anomaly in QCD and the chiral algebra: 

(tt+tt-I (£") 2 |0> = ^g 2 + 0{a s ql) + 0(m 2 ) , (7) 
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where q = p+ + P- is the total 4-momentum of the pion pair (so that qo is the total energy 
of the pair), b = 9 is the first coefficient in the QCD beta function with three light quarks, 
and the subleading terms, analyzed in Ref.^2], are relatively slowly varying with q 2 in the 
physical region of the pionic transition. In practice these terms can be approximated by a 
constant C, which depends on q Q , so that the amplitude of the decay can be written according 
to the equations (J2J), ©, and (|7J) as 

4vr 2 

A(2S - 7T+7T- IS) = -— a 1S -2S (q 2 ~ C) (8) 



with C being slightly different for the transitions in charmonium and in bottomonium. The 
observed experimental spectra of the dipion mass in the decays ip' — > tttt J/i/j and T' — > Tin T 
agree with the form of the amplitude in eq.(JE} and the with the numerical value of the 
constant: C = (4.6 ± 0.2) m 2 for the transition in charmonium and C = (3.3 ± 0.2) m 2 for 
the decay T' — > nn T (for a discussion see e.g. the review [T3]). 

The total rates of the pionic transitions are calculated from eq.(JH} as 

T(ijj' -> 7T+7T- J/ij>) w 0.30 \aj/4,-^\ 2 Ml/) - M(J/^)] 7 , 

r( T ' -> 7T+7T- T) « 0.36 |a T _ r ,| 2 [M(T') - M(T)] 7 , (9) 

where the decimal numerical factors describe the relative suppression of the rates due to the 
corresponding constant term C and due to the nonzero pion mass. Thus from a comparison 
of these expressions with the experimental data on the rates one estimates the transition 
chromo-polarizabilities: 

\aj/^ \ w 2.0 Gel/" 3 , |a T -T'| ~ 0.66 GeV^ . (10) 

These estimates of the transition matrix elements illustrate the typical values that one can 
expect for the diagonal chromo-polarizability in charmonium and bottomonium. It is also 
somewhat satisfying to notice that for the bottomonium the estimate of the diagonal chromo- 
polarizability from eq. (0) is in a reasonable agreement with the value of the transition term, 
since one generally would expect the diagonal matrix element © to be somewhat larger 
than the non-diagonal (JHJ)- 

Proceeding to discussion of the decays J/if> — > ttit£ + £~ and T — > tttt£ + £~ we retain the 
notation q for the total 4-momentum of the two pions. The amplitude of such decay for a l 3 ^ 
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state of heavy quarkonium in the soft pion limit can be written as a sum over intermediate 
n 3 Si states: 

a (1 > Si - ,wr) = \ (*>'!<» £ M(, ig )-M(i S) + ,/ <» 3s - - f+r > ■ 

(ii) 

where the sum goes over the discrete states as well as the continuum. In writing this 
expression it is taken into account that the soft pion approximation is only valid at go 
M(1S), so that any recoil of the heavy quarkonium upon emission of the pion pair can be 
and is neglected. In this limit the relation between g and the total momentum I of the 
lepton pair can also be written as M 2 (IS) - I 2 = 2 q M(1S). 

In the chiral limit the first term (with n — 1) in the sum in eq. ([TTf dominates for 
soft pions, due to its singular behavior as 1/qo- It can be noticed however that the decay 
amplitude itself is not singular due to the (even faster) vanishing of the pion production 
amplitude (J7J) in the limit of soft pions. In the 'real life' the minimal practical energy g 
is not much less than the spacing of the quarkonium levels, and the contribution of higher 
states mixes into the amplitude. In what follows we first consider the contribution of only 
the first term of the sum in eq.ljlljl and then discuss the effect of the higher terms. Keeping 
only the contribution of the first term in the sum in eq. (jll|) one can write the differential 
rate of the discussed decay in terms of the chromo-polarizability aig and the leptonic width 
r ee (l 3 51) = r(l 3 5! -> £+£') in the form 



[q 2 -C(q )] 2 ^ 2 / 4ml 
Ab 2 q 2 Ksl V ^ 



dY{l 6 S x -> n + n-£ + £-) = * \a is \ z \\l - V ql ~ q 2 T ee (l A Sl) dq 2 dq , (12) 



where the formula (|ZJ) is used with a simplified parametrization of the subleading terms as a 
constant C(qo) similar to that in eq.(jHJ- 

In order to assess the feasibility of observing the discussed decays it can be noted that 
at a given constraint on the maximal value of q$: q^ < A (or equivalently at a lower cutoff 
on the invariant mass of the lepton pair) the probability described by eq.(|12|) strongly peaks 
near the highest values of both q 2 and q$, i.e q 2 ~ A 2 and go ~ A, and the total probability in 
the kinematical region constrained as go < A scales approximately as A 6 . However at higher 
q 2 both the dominance of the diagonal IS — IS transition in the sum in eq. (jll|) becomes 
weaker and the linear in q 2 behavior of the amplitude in eq.([7|) derived for soft pions becomes 
questionable. It is still quite likely that with these limitations the presented here approach 
can be used up to somewhat higher values of A: A ps 0.8 — 0.9 GeV, than those observed 
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in the pionic transitions from if)' and T'. Indeed, experimentally the linearity in q 2 of the 
amplitude of the transition T' — > ixix T is very accurate |15t ITU] in the physical region, i.e. up 
to y/q 1 fa 0.56 GeV with no obvious hint of its violation close to this region. On the other 
hand the /o(980) resonance places a natural upper bound on the region of applicability of 
eq.(jZJ). As to the contribution of higher quarkonium states in the sum in eq. (|lljl . for each of 
these states the magnitude of this contribution relative to that of the diagonal transition is 
given by 

C^lS-nS <Zo 

Tn ~ a ls M{n 3 S 1 )-M{VS 1 ) + q 

The transition polarizability ais~ n s should considerably decrease with n. This is supported 
by the very small experimental rate of the decay T(3S) — > ixix T 1 . Thus, most likely, the 
only real effect of higher states up to A ~ 0.9 GeV reduces to that of the 2S resonances. 
These however can be accounted for in the data analysis, since for these resonances all the 
parameters (except for the overall relative phase of their contribution) in eq.(|13|) are known. 
Furthermore an observation and analysis of the discussed decays at higher values of q$, where 
the expression (jZJ) is no longer valid, would be of a great interest for studies of the pion-pion 
scattering beyond the soft-pion region. 

Numerically one can estimate from eq. (|12jl the total rate in the kinematical region con- 
strained by A = 0.9 GeV as 

2 r ee (l 3 Sa) . (14) 

Given that the diagonal polarizability is likely to be larger than the transition one, it can be 
expected that for the J / if) resonance the branching ratio of the discussed decay in a useable 
kinematical range should be at the level of 10~ 5 which looks to be well within the reach with 
the expected CLEO-c data sample. For the T resonances the effect is reduced by a factor of 
about 20 due to a smaller chromo-polarizability (cf. eq.QlUJl) and also due to a smaller value 
of B(T — > £ + £~). Thus an experimental study of the discussed decay for the T resonance is 
likely to be a future task for a 10 35 upgrade of KEKB[16J and/or the 10 36 B-factorv[T7j. 

lr Tlie known problem with this decay is that the observed dipion mass spectrum does not agree with the 
formula in eq.J7J) (see e.g. in the review A natural although quite qualitative explanation of this fact 

is that the polarizability a\s-3s is very small, and subleading effects in the multipole expansion in QCD 
come into play in this transition. 
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r ee (i 3 Si) 



(13) 



r(i 3 Si -> ix + ix- 



\q <0.9GeV 
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